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ABSTRACT 

Tilt angles of close to 30,600 sunspots are determined using Mount Wilson 
daily averaged magnetograms taken from 1974 to 2012, and MDI/SoHO magne- 
tograms taken from 1996 to 2010. Within a cycle, more than 90% of sunspots 
have a normal polarity alignment along the east-west direction following Hale's 
law. The median tilts increase with increasing latitude (Joy's law) at a rate of 
~ 0°.5 per degree of latitude. Tilt angles of spots appear largely invariant with 
respect to time at a given latitude, but they decrease by ~ 0°.9 per year on 
average, a trend which largely reflects Joy's law following the butterfly diagram. 
We find an asymmetry between the hemispheres in the mean tilt angles. On 
average, the tilts are greater in the southern than in the northern hemisphere for 
all latitude zones, and the differences increase with increasing latitude. 

Subject headings: Sun: dynamo - Sun: general - Sun: interior - Sun: magnetic 
fields - Sun: rotation - Sun: sunspots 

1. Introduction 

The sunspot magnetic tilt angle is the angle between the east-west direction and the 
line connecting opposing polarity portions in a sunspot group or in a pair of sunspots. It 
can be described by two well-known but qualitative laws: Hale's and Joy's laws. Hale's law 
states that, at any moment, sunspot polarity pairs have opposite sign orientation between 
the northern and southern hemispheres. At any given cycle, a majority of the leading spots 
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have opposite sign polarities in opposite hemispheres. The sign orientations switch at the 
onset of a new cycle when new sunspots emerge from high latitudes ( Hale fc Nicholson||1925 ). 
Joy's law describes the latitude dependence of the tilt of bipolar sunspot regions: the leading 
spots are closer to the equator than the following spots and this tilt angle increases with 



increasing latitude, as found in a study by Joy reported by Hale et al. (1919). One motivation 



of the present study is to obtain quantitative and accurate, long-term measurements of the 
tilt angle as functions of time and latitude. 

Accurate, long-term measurements of the tilt angle can be used to discriminate between 



published models. In the Babcock-Leighton model (Babcock 1961 Leighton 1969), which 



describes the repetitive solar global magnetic field patterns, the tilt angles of bipolar magnetic 
regions are interpreted as pitch angles of magnetic field lines stretched by differential rotation. 
This interpretation by itself does not reproduce Joy's law. On the other hand, stretching 
magnetic field by differential rotation in the latitudinal and radial directions is the basis of the 
magnetic buoyancy needed for flux tube emergence ( Parker|[l955a[b ) . Numerical simulations 
show that magnetic buoyancy and the Coriolis force together generate the correct latitudinal 
tilt angle distribution QD'Silva fc Choudhuri||l~993~l |Fan et d.]l994| |Fisher et al.|[l995| . These 
models also predict a close correlation between the magnetic field strength and the tilt angle. 



While some observations confirm the existence of such a correlation (Tian et al. 2003; Dasi- 



Espuig et al. 2010), others show no evidence for it. For example, a study of the tilt angles 



of over 700 sunspot groups during their disk passage found that the sunspot tilt angles are 
established at an early stage of the sunspot emergence and remain roughly constant until the 



decay of the active regions (Kosovichev & Stenflo 2008). A recent numerical simulation of 



a thin flux tube in a rotating turbulent spherical shell needs fewer constraints on the initial 



magnetic strengths and shows a statistical trend of the tilt angles with latitude (Weber et 



al. 2011). An initial flux strength near 40 kG produces tilt angles following Joy's law, while 



stronger and weaker flux tube strengths give a wider spread of tilt angles on the surface 
of the Sun. The differences between the Babcock-Leighton model and the magnetic flux 
buoyancy process underline the complex nature of the magnetic flux emergence as part of 
the solar dynamo operation. 

We measured the tilts of sunspots that erupted over the course of 38 years (1974-2012) 
using magnetic field data. Our analysis differs from an earlier one by Howard ( |1994 ) which 
was based on the white light images from Mt. Wilson from 1917 to 1985 from which sunspot 



positions were measured by Howard et al. (1984). In their study and others based on it 



magnetic polarities were not known and the position of the polarity reversal was assumed to 
be midway between the leading and trailing spots. The second difference between our work 
and earlier studies is that we examined the average tilt angles of a sunspot group using one 



magnetogram per day, while others studied daily sunspot tilt angle variations (Gilman & 
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Howard||1986 Howard] 199 laf b Muneer fc Singh]20 02). The time period that we covered was 



exactly following that period covered by Dasi-Espuig et al. (2010) for the tilt angle study. 



The various studies are also distinguished by their relative sizes: ours is larger than those 



studied by Wang & Sheeley (1989) and Kosovichev & Stenflo (2008) because our time period 



is longer, but smaller than the sample in a recent study by Stenflo & Kosovichev (2012). 



Their sample "spans a range of scale sizes across nearly four orders of magnitude" including 
sunspots. It is interesting that they found large-scale bipolar regions having tilt angles about 
90° from the average Joy's law distribution. These regions are larger than sunspot regions, 
are not included in our study. 

What mostly separates our study from many others using magnetograms is that we 
analyzed only identified sunspot and sunspot groups. By directly overlapping positions of 
sunspot groups on the magnetograms, we avoided the need to develop a more sophisticated 
algorithm for pattern recognition. We constructed magnetic bipolar regions by calculating 
centers of mass of the positive and negative polarities within a sunspot group. This approach 
results in the averaged tilt angles of a sunspot region, but not tilt angles due to local dis- 
turbances and fine magnetic structures. This is different from the tilt angles obtained over 
small-scale, quiet-Sun bipoles, or ephemeral regions comparable to the X-ray bright points 



(Sattarov et al. 2002 Tlatov et al. 2010) 



2. Data 

The data used in this study consist of sunspot daily records and longitudinal magne- 
tograms. For the period 1974 to 1990, we use the sunspot records made at MWO, which 
has assigned its own sunspot group numbers since 1920. The Mt. Wilson sunspot records 
do not contain the areas of the sunspot groups, but have the crucial information needed 
for this study, namely, the sunspot central meridian crossing time and latitude. From 1991 
onward, the sunspot reports were prepared by Space Weather Forecast Center (SWFC) from 
observations made by MWO and other observatories. They are made available via their 
web site. Parameters relevant to this study include the observing dates, sunspot NOAA 
numbers, sunspot group areas, and sunspot group disk locations. They correspond to the 
sunspot status valid at 00:00 UT each day. 

The magnetograms for the entire study period were obtained at the 150 foot solar tower 
telescope at Mount Wilson Observatory. A two-mirror coelostat is installed on top of the 
tower, reflecting light to the objective lens with a focal length of 150 feet. Solar images are 
formed at ground level in an observing room. A spectrograph is situated under the ground 
level in a vertical pit. Sunlight enters the spectrograph through a slit, and exits into one 
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of several fiber-optic bundles at selected wavelengths. The first Mt. Wilson magnetograph 



was built by Babcock (1953) and daily magnetograms of the full disk of the sun were started 



in 1957. The magnetograph was largely rebuilt by Howard et al. (1983). It consists of 



an interference filter having a 100 A bandwidth centered on the Fe I A5250 line, a KD*P 
crystal-Glan Thompson prism combination providing circular polarization modulation; a 
Littrow spectrograph with a new grating; and a new exit slit assembly. The current system 



was an upgrade to a 24-channel spectrophotometer by Ulrich et al. (1991, 2002) adding 



three spectral lines: Cr II A5237, NaD A5896 and Nil A6768. Since the mid-1980s, daily 
observations at MWO generate the magnetograms and intensity-grams taken in the 5250.2 
A Fe I and 5237.3 A Cr II lines. These data are used to model the total solar irradiance 



variations (Ulrich et al. 2010). 



The systematic upgrade of the magnetograph at MWO, while largely preserving the 
spectral sampling characteristics, provides consistent magnetic field observations over many 
decades. The full-disk daily-averaged magnetograms at the spectral line A5250.2 A have the 
longest observing history, and are used in this study. A daily magnetogram was constructed 
by averaging as many as 20 images corrected to the observing time at 20:00 UT. Differential 
rotation is taken into account when the images are combined. The non-uniform pixel size 
due to the nonuniform image sampling during the solar image scan was corrected. The 
final image size is 512 x 512 pixels, and the pixel size is 3.7" x 3.7". The daily averaged 
magnetograms are available from 1985 to the present time. Prior to 1985, the individual 
magnetograms obtained once a day were used for the study. 



From 1996 January to 2010 December, we employed the full disk MDI/SoHO (Scherrer 



et al. 1995) magnetograms, level 1.8, with a 96-minute time cadence. The MDI employs a 



Michelson interferometer as the tunable spectral device, centered at Ni I 6768 A. The image 
size is 1024 x 1024 pixels, and the pixel size is ~ 2" x 2". A single MDI magnetogram per 
day was employed for all sunspots appearing on that day for the tilt angle calculations. That 
magnetogram was taken closest to 00:00 UT when the sunspots were recorded. 



Method 



Two definitions representing different aspects of the tilt angle are useful for studying 
Hale's and Joy's laws. They are illustrated in Fig. 0. Fig. ([IJi) is similar to that defined 



by Howard (1991b). The tilt angles, 7, vary from < 90° to > 90° as the positive polarity 
component changes from the leading to the trailing component. This definition carries the 
information of the leading spot polarity signs. Hale's law is conveniently illustrated using 
this definition. The orientation of tilt angles in Fig. ([ljo) is identical to those defined in 
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Fig. but the bipolar magnetic polarities are no longer reflected by the values of tilt 



angles. The tilt angles vary between [—90°, 90°] (Wang & Sheeley 1991), and are the most 
suitable to illustrate Joy ? s law. It is noted that signs of the tilt angles are the opposite of the 
conventional angles in a Cartesian coordinate system in Fig. (fib) . This ensures that Joy's 



law will be conveniently described as the tilt angles increasing with latitude. 

The tilt angle measurements are performed by an automated IDL package which ac- 
complishes the following tasks: (a) locate the sunspots on the magnetograms; (b) determine 
the centroids of the entire sunspot, positive and negative polarities within each spot; and (c) 
calculate the tilt angles for each sunspot group. Once the sunspot locations and sizes were 
initially determined, the automated IDL program runs an iteration of tasks (b) and (c) until 
the tilt angle converges or stabilizes. 

Specifically, we first employed an ellipse to define the boundary of an active region. 
The center of an elliptical area is the sunspot location tabulated by MWO and SWFC. 



The disk locations of the sunspots are corrected for differential rotation (Newton & Nunn 



1951) in order to match the sunspots with the magnetograms. The orientation of the ellipse 
representing the tilt angle is set to 0° because most sunspot groups are elongated in the E-W 
direction. After some experiments, we found that reasonable tilt angles were produced if the 
long and short semi-axes, a and b, of the initial ellipse were set as 



a = 20 x (sunspot area) ' , b = a/2 (1) 

where sunspot area is the area given by SWFC in millionths of the solar hemisphere, but was 
converted to image pixels. For sunspots prior to 1991, sunspot area is set to 60 millionth 
hemisphere. We note that the value of the initial area does not appear critical for the tilt 
angle calculations as long as the sunspot locations are roughly correct. At least, this is 
true to the MWO magnetograms. Empirically, the search program converges to the region 
surrounding the sunspot groups. 

Following the initial setup of the ellipse over a sunspot group, a new centroid of the 
active region is computed with the mass center of the total magnetic flux and the pixels are 
sorted according to their polarity. From these two pixel sets, we find a pair of circles defining 
a simple bipolar configuration, which radii are defined as \J area/ir, where area is the total 
number of pixels in each polarity. The idealized bipolar regions often have overlapping circles 
of opposite polarity because pixels of either strong or weak magnetic fields are given equal 
counts to sum up the areas. In reality, the strong opposite polarities appear side by side in 
an active region, but weaker field spreads throughout the region. In our analysis, magnetic 
field signals are extracted down to the minimum field strength of 10 G for MWO data, and 
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20 G for MDI data. Figure ([2]) shows examples of three cases of separated, contacted and 
overlapping magnetic polarity regions fitted by our algorithm. 

The measured centroids of the positive and negative polarities are in the observed solar 
disk coordinate system designated as [x", y"}. We convert them into heliographic latitude and 
longitude, [B, L] degrees, taking into account the instantaneous heliographic latitude, B . 
The conversion is available with the code xy2lonlat which is implemented in the SolarSoftware 

by Thomas Metcalf. An updated magnetic tilt angle is calculated 
from 7 = arctan(A5, AL x cos 5), where AB and AL are the differences between positive 
and negative polarity centroids in heliographic latitude and longitude; B is the heliographic 
latitude of the positive polarity centroid. 

The iteration continues with a new set of parameters: a non-zero tilt angle, a centroid of 
the total magnetic flux, and centroids of positive and negative polarities; the long semi-major 
axis was calculated by newly determined parameters, 

a = 3(r p + r n + s) (2) 

where r p and r n are the radii of the circular areas of idealized positive and negative polarities; 
and s is the separation of the bipole. The semi-minor axis, b, is the bigger of values between 
r p and r n . Successive tilt angles were compared, and the iteration is set to stop when the 
difference is less than or equal to 0°.l. Convergence takes <20 iterations for 70% of sunspots 
while ~80% of spots converge to 0°.l within 24 iterations. Visual inspection confirmed that 
the tilt angle determinations are reasonable for even the largest number of iterations. The 
process is illustrated in Fig. ^ where multiple ellipses overlapping a sunspot group on a 
magnetogram show successive solutions from the algorithm. The initial ellipses are evident 
by their zero orientation with respect to the equator. The fitted ellipses quickly turn to 
the final orientations during the process (black ellipses). The final derived tilt angles of the 
sunspots in Fig ^ are shown in Fig. Q. 

The fact that we start with identified sunspots in step (a) minimizes confusion in the 
identification of active regions on the magnetograms. Without anchoring the method with 
pattern recognition that originally came from human observers, the algorithm would have 
required a much more sophisticated method. 

In our analysis, only sunspots having central meridional angle < 45° were used for the 
tilt angle study in order to minimize errors caused by projection. This is equivalent to a 
6-day period of a distinct sunspot group disk passage. The tilt angles were calculated with 
the automated IDL program for about 28642 sunspot groups observed by MWO from 1974 
to 2012 March and 11932 by MDI from 1996 to 2011 April. Although the MDI data have a 
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shorter observing time span (1996-2011) compared to MWO data, MDI has higher observing 
rates than MWO because of the unlimited weather conditions and non-local time zone. 

About 20% of sunspots in MWO data, and 12% of sunspots in MDI data were considered 
to have poorly determined tilt angles, and were visually re-examined after the automated 
program run. These sunspots were singled out by the following criteria: (a) During the 
disk passage, those sunspot groups had standard deviations of tilt angles > 20°. 5. This is 
equivalent to a tilt angle varying from 60° to 120° for a sunspot group in 6 days. We note that 
many sunspot groups have tilt angle variations greater than 20°. 5 in the course of maximum 
6 days, but yet the tilt angles appear valid. Therefore, 20°. 5 is a low threshold to detect the 
ill-defined tilt angles; (b) Those sunspots having latitudes less than 15° but the tilt angles 
greater than 45°. This ensures that a large tilt angle is valid for sunspots at low latitudes; 
(c) Those sunspots having bipolar separation greater than 15°. Based on the SWFC sunspot 
reports from 1991-2012, about only 1.5% of sunspot groups have longitudinal extent greater 
than 15°. 

For those sunspot groups with ill-defined tilt angles, we re-determined the tilt angles 
by manually entering the sunspot central positions, the semi-major and minor axes and the 
initial tilt angles. During the manual process, some sunspots were rejected altogether due to 
the poor observations of sunspots on a magnetogram or they were too close to other major 
sunspots. About 2000 sunspots were rejected from MWO data, and about ~ 150 sunspots 
were rejected from MDI data. These correspond to 7.5% of MWO and 1.2% of MDI data 
samples, respectively. Although we rejected these sunspots, the following results would not 
be materially altered by their inclusion. 

The MWO and MDI data overlapped from 1996 to 2010, allowing us to compare results 
from these independent data sets. We find that the sunspot magnetic tilt angle measure- 
ments agree remarkably well between two different instruments. This is shown in Fig. ^ 
where tilt angles from the two instruments are distinguished by colors, red (MW) and blue 
(MDI). The sunspot butterfly diagram is shown in the background in which latitudes are 
centroids of sunspot groups measured from MWO (dark colored "+" ) and MDI (light colored 
"+"). Because of the agreements between tilt angle measurements from two instruments, we 
combine them for the time period from 1974 to 2012. Specifically, tilt angles presented in 
the following sections were measured from MWO data between 1974 to 1995, and after 2010. 
The tilt angles were measured from MDI between 1996 and 2010. MDI data gaps between 
June 1998 to January 1999 were filled with the tilt angle measurements from MWO. Over 
all, the total sunspot number used for the tilt angle study is 30623, and they made up 8705 
distinct sunspot groups. 
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4. Results 
4.1. Hale's Law 

The tilt angles of all 30623 sunspots are presented in Fig. ([6]). Hale's law is evident 
in that the majority of sunspots have opposite polarity distributions in the northern and 
southern hemispheres at any given cycle. The leading polarity switches signs in the period 
when an old cycle fades, and a new cycle rises. We tabulate the percentages of sunspot 
numbers with their leading polarities either positive (7 < 90°) or negative (7 > 90°) in the 
respective hemisphere from cycle to cycle. We adopt the official cycle definition available at 
the FTP sit<Q In our analysis, starts of cycles are those in the column "Year of Minimum" , 
and the lengths of cycles are those listed in the column "Cycle Length". Three complete 
cycles were covered by our study, 21, 22, and 23. Within each cycle, Table ([!]) shows that 
more than 90% of sunspots had "normal" polarity alignment along the east-west directions. 



4.2. Joy's Law 

Fig. Q shows histograms of the sunspot tilt angle distributions in four latitudinal 
zones in the northern (red) and southern (blue) hemispheres. The tilt angles were measured 
between [—90°, 90°] using the definition illustrated in Fig. (JTJd) . Within each latitude zone, 
the sunspot numbers peak at different tilt angles in the two hemispheres. The separation 
between the twin peaks increases with increasing latitude. 

The median tilt angle, is the angle dividing the sunspot tilt angle distribution function, 
n s (7), into equal areas. To compute 7, we construct the tilt angle cumulative distribution 
function from n s {y), 

N,(y) = n s ( 7 )d 7 (3) 

J -90° 

The median tilt angle is half the maximum cumulative function, 

7 = F(imax[JV a ( 7 )] > ) (4) 



2 

where F is the inverse function of N s (y). The uncertainties of the median tilt angles were 



1 ftp: / / ftp.ngdc.noaa.gov/STP /SOLARJDATA/SUNSPOT JSfUMBERS /INTERNATIONAL / maxmin/M AXMIN 
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estimated from 3 x ^/\fW~ s , where N s is the total number of sunspots used to calculate 
the median tilt angle. Table (|2j) lists the median tilt angles calculated from the tilt angle 
distributions. It is consistent with what is shown in Fig.(j7j). The median tilt angles clearly 
increase in absolute value with increasing latitude. 

A tilt angle hemispheric asymmetry is shown in Fig. ([7]) and Table (|2]). The sunspot 
number peaks closer to 0° tilt angle in the northern hemisphere than in the southern hemi- 
sphere. We found that the asymmetry in the tilt angles between the hemispheres existed 
in all cycles. On average, the tilt angles were higher in the southern than in the northern 
hemisphere by ~ 5° in the low latitude zone (0° — 10°), ~ 10° in the intermediate latitude 
zone (10° - 30°) and ~ 13° in the high latitude zone (> 30°). 

Joy's law is also illustrated in Fig. (p]) in which median tilt angles were calculated using 
Equations ^ and Q within 5°-latitudinal bins. The linear correlation between the sunspot 
latitude, B, and the median tilt angle is written j(B) = kB + c. Both constants k and c are 
to be determined by fitting data points (black circles in Fig (|8|). The linear least-square fit 
produces the red straight line, 



7(5) = (0.5 ±0.2)5 - (0°.9±0°.3). (5) 

The uncertainty of the slope was propagated from the median tilt angle uncertainties. Two 
red straight dotted lines are fits to the median tilt angles ± error bars. The differences 
between two fits gave the errors of 3cr, which are ±0.2 for Joy's law slope, and ±0°.3 for the 
tilt angle at B — 0°. The tilt angle hemispheric asymmetry discussed above is the reason 
why the constant term, c = — 0°.9, but not zero. 



Wang & Sheeley (1991) gave an empirical equation for tilt angles by assuming that the 
Coriolis force acts on an emerging flux rope as it undergoes expansion in the longitudinal 
direction while it is below the sun's surface. Their study was based on examining 2700 



bipolar magnetic regions (BMRs) erupted between 1976-1986 (Wang & Sheeley 1989). The 
tilt angle, 7, of the flux tube follows sin 7 ~ u(B)r ex sin 5, where oj(B) is the differential 
rotation rate as a function of latitude B, r ex represents the sub-surface flux rope expansion 
time. For a typical solar rotation rate 14°.4/day, and r ex = 2.0 days, we obtain u(B)r ex = 0.5 
which is the correlation coefficient between sin 7 and sin B. Coincidently, this agrees very 
well with our measurements: sin 7 = (0.5 ± 0.2) sini? when sin 7 and sinB are fitted with a 
linear least square fit. 
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4.3. Tilt Angles and Sunspot Butterfly Diagram 

The average latitude of sunspots depends on the phase of the sunspot cycle as illustrated 
by the butterfly diagram. The tilt angles may also vary with time due to Joy's law. In 
order to determine if the time dependence of the tilt angle is strictly a consequence of 
the sunspot changing latitudes, we calculated median tilt angles, 7, using the histogram- 
cumulative function technique described by Equations ^ and Q. The sunspot tilt angle 
distribution function is constructed in a series of time intervals, t — At/2 and t + At/2, where 
At = 150 days in our analysis. The measurement uncertainties were calculated, again, as 3a, 
where a = ^f/\/W s . Fig. ([9J shows the median tilt angles as a function of time overlapping 
a butterfly diagram made from the centroids of the sunspot latitudes determined from the 
total magnetic flux. The latitude scale for sunspots is on the left vertical axis, and the tilt 
angle scale is on the right vertical axis. The median tilt angles are plotted as red curves, and 
circles. It appears that tilt angles decrease with time in both hemispheres as the sunspots 
migrate from high to low latitudes through each activity cycle. 



Fig (10) shows median tilt angles over the monthly averaged sunspot number (available 
at |http : / / solarscience . msf c . nasa . gov / greenwch / spot _num . txt \ . The tilt angles were averaged 
between two hemispheres shown in Fig ^ in order to increase the signal-to-noise. The 
uncertainties were calculated accordingly. The fact that tilt angles show the general declining 
trend with cycle phases allows us to combine the tilt angle measurements of cycles 21, 22 
and 23. We registered tilt angle measurements at times of cycle minima with equal time 
intervals, 3652 days (~ 10 years). This is illustrated by three horizontal bars on top of figure 



(10). With improved signal-to-noise levels, Fig (11) shows the combined tilt angles as a 
function of time within a typical cycle. Error bars represent ±1 standard deviations of tilt 
angles of three cycles. Large error bars of tilt angles from year to 1 were caused by fewer 
sunspots and their diverse tilt angles during the cycle exchange period. By fitting tilt angles 
(red circles) with an error- weighted linear least-square fit, we obtain 



j= (-0°.9±0°.l) t+ (11°.2±2°.0) (6) 

where t (year) is the time since the start of the solar cycle. The errors of averaged tilt angles 
produced the two red dotted lines. The differences between the lines are the uncertainties in 
Equation ([6]). On average, the tilt angles decrease ~ 0°.9 per year. Within a typical 11-year 
solar cycle, the tilt angles decrease about 10°. 



The latitudes of sunspot centroids averaged from three complete cycles are shown in 
Fig (12). It is consistent with the butterfly diagram that the sunspots generally migrate 
from high to low latitudes as the cycle progresses. The median sunspot latitudes and their 
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uncertainties as a function of time were linearly fitted: B(t) = (— l o .7±0°.2)t + (23 o .3±2°.7), 
where t is measured in year from the beginning of the cycle. On average, sunspot groups 
migrate from high to low latitudes at the rate ~ 1°.7 yr _1 , while the tilt angle decreases 0°.5 
per degree of latitude (see Equation ([5])). This results in the tilt angle decreasing roughly 
0°.9 yr _1 which rate is consistent with Equation This implies that sunspot tilt angles 
decrease with time following Joy's law. 



4.4. Tilt Angles and Cycle Phases 

Within a limited latitude range, the sunspot tilt angles were approximately invariant 



with respect to time. Fig. (13) shows the tilt angles averaged over three cycles within three 
latitude zones as functions of time. Spots with \B\ > 30° are not plotted because there are 
too few of them, and the uncertainties are too great. The data points are represented by 
colored circles and solid curves. Joy's law is, again, evident in that the median tilt angles 
are generally higher as the latitude zones increase. The blue curve connecting blue circles 
represent the sunspots at latitudes 20° < \B\ < 30° and stays generally at the highest 
tilt angle range. The green curve connecting green circles represent sunspots at latitudes 
10 < \B\ < 20° and occupies the intermediate tilt angle range. The red curve connecting 
red circles represent sunspots at latitudes \B\ < 10° and outlines the low tilt angles. Within 
the tilt angle uncertainties (not shown in the figure), the figure shows that tilt angles stayed 
roughly constant through time at a fixed latitude zone. We plot the median tilt angles with 
colored horizontal lines corresponding to each latitudinal zone. Above respective lines, the 
median tilt angles and their uncertainties are marked with respective colors. The sunspot 
tilt angles are independent of cycle phases. 

The seeming contradiction between observations in Fig. (11) or Equation ^ and Fig. 



(13) can be understood by numbers of sunspots varying with cycle phases and latitudes. 



In Fig. (13), three color-dotted curves represent sunspot numbers used to calculate the tilt 
angles at three different latitudinal zones (the scale is shown by the right vertical axis). On 
average, there were fewer sunspots at low latitude (red and green dotted curves) than at high 
latitude (blue dotted blue curve) in the beginning of the cycle (around year=l). This brings 
the median tilt angles calculated from Equations ^ and ^ to higher values. Likewise, 
the red and green dotted curves rise and the blue dotted curve declines with time. As the 
result, the median tilt angles decrease as sunspots migrate from high to low latitude as cycle 
progresses, but the median tilt angles are independent of cycle phases at a fix latitudinal 
range. 
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5. Summary 

We calculated the sunspot magnetic tilt angles using daily average magnetograms by 
the Mt. Wilson Observatory from 1974 to 2012, and MDI magnetorams from 1996 to 2010. 
The tilt angles were measured from 30650 sunspots. We summarize the statistical results of 
the sunspot tilt angles in the time period 38 years: 

1. Hale's law is evident in the measured sunspot magnetic tilt angles. Within each cycle, 
over 90% of sunspot polarity follow Hale's law in both hemispheres, and sign alignment 
is reversed by the start of each new cycle. 

2. On average, Joy's law is evident in that the tilt angles increase with increasing latitude 
following the relation 'y(B) = 0.5B — 0°.9, where B° is the latitude. 

3. On average, the tilt angles decrease with time at the rate 0°.9 per year. This largely 
reflects Joy's law following the sunspot butterfly diagram. 

4. Within latitudinal zones, the tilt angles were independent of the cycle phases. The 
average tilt angles are 2° in the low latitudinal zone (\B\ < 10°), 6° in the median- 
low latitudinal zone (10° < \B\ < 20°) and 12° in the median-high latitudinal zone 
(\B\ > 30°). 

5. We found an unexplained, persistent asymmetry between the median tilt angles mea- 
sured in the north and south hemispheres in all latitudinal ranges. The average tilt 
angles were greater in the southern hemisphere than in the northern hemisphere by 
~6° in the to 10° latitude zone, by ~10° in the intermediate latitudel zone (10° — 30°) 
and by ~13° in the high latitude zone (> 30°). 
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Table 1. Hale's Law: Leading spots with positive magnetic polarities 



Cycle (Length) 


Year of 


Northern Hemisphere 


Southern Hemisphere 


[yrs] 


Minimum 


SSN 


7 < 90° 


SSN 


7 < 90° 


21 (10.3) 


1976.5 


4425 


91.5% 


4655 


8.7% 


22 (10.0) 


1986.8 


3781 


9.1% 


4130 


91.5% 


23 (12.2) 


1996.9 


5457 


92.3% 


6272 


6.5% 



Note. - - "SSN" represents the total sunspot number used in the calculation 
of tilt angles in the respective categories. 7 is the tilt angle. Based on the tilt 
angle definition illustrated in Fig. ([T^.), the leading spots have positive magnetic 
polarities when 7 < 90°. 
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Table 2. Latitude Dependence of Tilt Angles (1974-2012) 



Latitude Northern Hemisphere Southern Hemisphere 
Zone SSN 7 ± 3a SSN 7 ± 3a 



< 10° 3449 

10° - 20° 7442 

20° - 30° 3374 

> 30° 519 



-2°.2±0°.l 3669 

2°.l±0°.l 7572 

7°.0 ± 0°.4 3888 

7°.0±0°.9 710 



-8°.3±0°.4 
-12°.2±0°.4 
-16°.5±0°.8 
-19°.9±2°.2 



N 

90* 



© 



E m ' (?) (S) @ 0' W 



® 

-90* 

(a) 



Tilt Angle >0 
trailing polarity 




Tilt Angle <0 

leading polarity 



leading polarity 




trailing polarity 



equator 



(b) 



Fig. 1— Definitions of the tilt angles, (a) Tilt angles vary between [-180°, 180°]. "P" 
represents the positive polarity, and "N" represents the negative polarity, (b) Tilt angles 
vary between [—90°, 90°]. The orientations of the tilt angles are identical in both definitions. 
The difference is that the leading spot polarity is indicated in the definition (a), but not in 
the (b). 



Fig. 2. — The ellipses represent the sunspot group boundaries. The bipolar spot pair is 
represented by red and yellow circles for positive and negative polarities. Three extreme 
cases are illustrated with polarity pair circles separate, contacted and overlapping. The long 
axes of the ellipses with respect to the equator represent the magnetic tilt angles of each 
bipole. Magnetic field strength is displayed on a scale from -250 (black) to 250 G (white). 
The size of each panel is 320x200 pixels (635" x 398"). East is to the left, and the north is 
up. The magnetograms are from MDI data. 
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Fig. 3. — Illustration of the tilt angle measurement algorithm. The multiple ellipses around 
each sunspot represent successive iterations of the fitting algorithm. The initial ellipses are 
parallel to the equator. They are quickly adjusted to the orientation close to the final state, 
shown as black ellipses. 
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Fig. 4. — Measured tilt angles are represented by the tilts of the long axes of ellipses with 
respect to the equator. Two circles inside the ellipses represent the mass centers and sizes of 
positive (red) and negative (yellow) magnetic polarities of sunspots. Equator is represented 
by the thick arc curve marked EW. These are the same regions as those in Fig. (J3T) 
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1996.0 1997.4 1998.8 2000.2 2001.6 2003.0 2004.4 2005.8 2007.2 2008.6 2010.0 

Year 

Fig. 5. — Sunspot magnetic tilt angles were measured from MWO (red curve) and MDI 
(blue curve) magnetograms. They are over-plotted with latitudes of sunspot mass centers. 
The dark colored "+" represent sunspots measured with MWO, and the light colored "+" 
represent those with MDI. The horizontal line represents the equator and 0° tilt angle. The 
scales for the sunspot latitudes is to the left vertical axis, and scales for the tilt angles is to 
the right vertical axis. The horizontal axis is time measured in "year" . 
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Fig. 6. — Hale's law during 1974 to 2012. The time spans from the end of cycle 20 to the 
beginning of cycle 24, and is marked in year. The horizontal line represents the equator. 
Red dots represent sunspots in the northern hemisphere, and blue dots represent those in 
the southern hemisphere. 
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Fig. 7. — Sunspot tilt angle distribution in northern (red lines) and southern (blue lines) 
hemispheres within four latitude ranges. The histogram bin size is 5°. The vertical straight 
lines mark the 0° tilt angle (see Table [2]) . 
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Fig. 8. — Illustration of Joy's law. Circles and solid curve represent median tilt angles 
of all 30,600 sunspots calculated with Equations ^ and Q in the 5° latitude bins. The 
vertical bars represent uncertainties calculated with 3 x j/^ssn, where ssn is the total 
sunspot number for the median 7. The solid and dotted red lines are the linear least square 
fits to the data points and their uncertainties respectively. Joy's law is expressed as 7 = 
(0.5 ± 0.2)5 - (0.9 ± 0.3), where B is the latitude. 
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1974 1977 1980 1983 1986 1989 1992 1995 1998 2001 2004 2007 2010 2013 

Year 

Fig. 9. — Median tilt angles (red curves) varying with time [year] overlapping the butterfly 
diagram. The median tilt angles were calculated with Equations (|3| and Q in the 300 
day interval. The uncertainties were estimated as 3 x 'j/^ssn, where ssn is the number of 
sunspots which contributed to the median tilt angle. The butterfly diagram is made with 
centroids of sunspots. The latitudes are scaled in the left vertical axis, and the tilt angles 
are scaled in the right vertical axis. The horizontal straight line indicate the latitude and 
tilt angle 0°. 
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1974 1977 1980 1983 1986 1989 1992 1995 1998 2001 2004 2007 2010 2013 



Year 

Fig. 10. — Median magnetic tilt angles (red curve and circles) and sunspot monthly average 
numbers (black curve) as functions of time [year]. The tilt angles were combined from both 
hemispheres (shown in Fig. [9]), and uncertainties were calculated accordingly. The scale for 
tilt angles is to the left vertical axis, and the scale for the sunspot numbers is to the right 
vertical axis. Three horizontal bars on top of the plot represent the starts and durations of 
cycles 21, 22, and 23. The starting dates and durations are listed in Table Q. 
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Year 

Fig. 11. — The median tilt angles (red) and average sunspot numbers (black) as functions 
of time (year) within a typical cycle. The tilts and sunspot numbers are averaged over 
three cycles, 21, 22, and 23 which were registered at the respective cycle minima with equal 
time intervals 3652 days ( ~ 10 yrs). The uncertainties of the tilt angles were the standard 
deviation of tilt angles at each data points. The tilt angles are fitted with the error-weight 
linear least square fit (solid red line) 7 = (— 0°.9±0°.l)t + (11°.2±2°.0), where t is measured 
in year from the beginning of the solar cycle. Two red dotted lines represent the uncertainties 
of fittings. 
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Fig. 12. — The sunspot centroid latitude (red) and average sunspot numbers (black) as 
functions of time (year) within a typical cycle. The tilts and sunspot numbers are averaged 
over three cycles, 21, 22, and 23 which is the processed in the same way as those in Fig. 
(11). The uncertainties were 1.5 times standard deviations of latitudes at each data points. 
The linear least square error- weighted fit is plotted by the solid red line, and two dotted 
lines represent the uncertainties of fittings. The relation is given by B = (— 1° .7 ± 0°.2)t + 
(23°. 3 ± 2°. 7), where t is measured in year from the beginning of the solar cycle. 
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Fig. 13. — Median tilt angles (solid curves, and circles) and average sunspot numbers (dotted 
curves) within three latitude zones as functions of time. The scale for the tilt angles is to the 
left vertical axis, and the scale for the sunspot numbers is to the right vertical axis. The tilts 
are averaged over three cycles, 21, 22, and 23. Three average tilt angles (colored horizontal 
lines) and their standard deviation at the respective latitude ranges are written above the 
colored lines. Colors represents latitude zones. Red: \B\ < 10°; Green: 10° < \B\ < 20°; 
and Blue: 20° < \B\ < 30°. 



